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1. Introduction 



Theoretical and experimental studies of the heavy quark sector in the Standard Model are 
of a great interest to complete the whole quantitative picture of fundamental interactions. 

In the bottom quark physics experimentalists step from the 10 6 yield of hadrons containing 
6-quark at present facilities up to the 10 9 yield in the foreseeable future to measure rare 
processes like the CP-violation and possible effects beyond the SM. To distinguish the hadronic 
dynamics from the latter effects at the quark level one needs the perfect understanding of QCD 
interactions binding the quarks into hadrons. 

One of the accompanying problem is the observation and study of the (be) state, yielding 
a 10~ 3 fraction of beauty hadrons at high energies. 

The B c meson allows one 

1. to accomplish the QCD-based models of hadrons with the bottom quarks, 

2. to study the specific production and decay mechanisms, 

3. to measure the SM parameters. 

The basic state of B c is the long-lived heavy quarkonium, which can be searched for in a 
way analogous to the observation of beauty mesons with a light quark [[J. 
At CDF the background is still strong to isolate the B c event at low statistics available ||. 
At present, the LEP Collaborations have reported on several candidates for the B c decays ||. 

Table 1: The basic characteristics of events with the 5 c -meson candidates at LEP. 



Collab. 


Mode 


mass of B c , GeV 


lifetime, ps 


OPAL,*, 




6.31 ±0.17 




OPAL new) 


lf)7l^~ 


6.29 ±0.17 


-0.06 ±0.19 


OPAL new) 


lf)lT^~ 


6.33 ±0.06 


0.09 ±0.10 


ALEPH 


1)1*1/ 


5.96 ±0.24 


1.77 ±0.17 


DELPHI 




6.35 ±0.09 


0.38 ±0.06 


DELPHI 


^(3tt) ± 


6.12 ±0.02 


0.41 ±0.07 



The mean values averaged over the i/jtt mode are equal to 

m Bc = 6.33 ±0.05 GeV, (1) 
r Bc = 0.28i°i° ps. (2) 

OPAL has reported also 

f{b -> B+) ■ BR(B+ -> i)7i + ) = (3.8t|£ ± 0.5) • 10" 5 . 

The question of the report is 

What are the theoretical expectations of the B c -meson properties? 
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2. The mass spectrum of the (be) family 



The most accurate estimates of (be) masses [HI [5] can be obtained in the framework of 
nonrelativistic potential models based on the NRQCD expansion over both I/uiq and v re i — > 

The uncertainty of evaluation is about 30 MeV. The reason is the following. The potential 
models |7J were justified for the well measured masses of charmonium and bottomonium. 
So, the potentials with various global behaviour, i.e. with the different r — > oo and r — > 
asymptotics, have the same form in the range of mean distances between the quarks in the 
heavy quarkonia at 0.2 < r < 1 fm ||. The observed regularity is the distances between the 
excitation levels are approximately flavor- independent. The latter is exact for the logarithmic 
potential (the Feynman-Hell-Mann theorem), where the average kinetic energy of quarks T 
is a constant value independent of the excitation numbers (the virial theorem) ||. A slow 
dependence of the level distances on the reduced mass can be taken into account by the use of 
the Martin potential (power law: V(r) = A(r/r ) a +C, a < 1) [|IIJ], wherein the predictions 
are in agreement with the QCD-motivated Buchmuller-Tye potential with the account for the 
two-loop evolution of the coupling constant at short distances [[11 . 

So, one gets the picture of (be) levels which is very close to the texture of charmonium and 
bottomonium. The difference is the jj-binding instead of the LS one. 

The spin-dependent perturbation of the potential includes the contribution of the effective 



one-gluon exchange (the vector part) as well as the scalar confining term |12 



Tr ,^ (L-S c L-S b \ ( dV(r) 8 1 \ 

+ - a s ^-+o a so b c - b b 4ir b{r) (3) 

3 m c m b r 6 3 6m c m b 

41-- - - 1 r 

+- a s (3(£ c • n) (S b ■ ft) - S c ■ S b ) — , n = - . 

3 m c m b r A r 

The model-dependent value of effective a s [Q can be extracted from the data on the 
splitting in the charmonium 

M(V>) - M( Vc ) = a s -^\R(0)\ 2 w 117 MeV. 

We take into account the renormalization-group dependence of a s at the one-loop accuracy 
by means of introduction of the quarkonium scale j| 

V 2 = (P 2 ) = 2(T)m red . 

The estimated difference between the masses of basic pseudoscalar state and its vector exci- 
tation [[| is equal to 

M(B* + ) - M(B+) = 65 ± 15 MeV. 
The mass of the ground state J| equals 

M(B+) = 6.25 ± 0.03 GeV, (4) 

which is inside the quoted region of the B c candidates at LEP. 
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Figure 1: The mass spectrum of (be) with account for the spin-dependent splittings. 
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Table 2: The masses of bound (fec)-states below the threshold of decay into the pair of heavy 
mesons BD, in GeV, in the models with the Martin and BT potentials. The spectroscopic 
notations of states are n 2jc Lj, where j c is the total angular momentum of c-quark, n is the 
principal quantum number, L is the orbital angular momentum, J is the total spin of meson. 



state 


Martin 


BT 


I 1 S 


6.253 


6.264 


l 1 ^ 


6.317 


6.337 


2% 


6.867 


6.856 


2 1 S 1 


6.902 


6.899 


2 1 P„ 


6.683 


6.700 


2P 1+ 


6.717 


6.730 


2P 1'+ 


6.729 


6.736 


2 3 P 2 


6.743 


6.747 


3^0 


7.088 


7.108 


3P 1+ 


7.113 


7.135 


3P V + 


7.124 


7.142 


3 3 P 2 


7.134 


7.153 


3D 2~ 


7.001 


7.009 


3 5 P> 3 


7.007 


7.005 


3 3 Pi 


7.008 


7.012 


3P 2 /_ 


7.016 


7.012 
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Radiative transitions 

The bright feature of the (be) family is that there are no annihilation decay modes due to the 
strong interaction. So, the excitations, in a cascade way, decay into the ground state with the 
emission of photons and pion-pion pairs. 

The formulae for the El-transitions are slightly modified. 

T(nP J ^n 1 S 1 + 1 ) = ^cwQeW I 2 (nP;nS) Wj(nP) , 

T(nP J ^n 1 S + 1 ) = ^-a em Q 2 cfi uj 3 I 2 (nP;nS) (1 - wj(nP)) , 

T(n 1 S l ^nP J + 1 ) = ^a cm Q 2 cS u 3 I 2 (nS;nP) (2J + 1) Wj (nP) , (5) 

T(n 1 S ^nP J + 1 ) = ^-a em Q 2 efi u 3 I 2 (nS;nP)(2J+l)(l-wj(nP)), 

y 

T{nPj -> nD, r + 7 ) = — a em Q 2 cS u 3 I 2 (nD; nP) (2 J' + 1) 

Wj(nP))wjr(nD)Sjj> , 
T(nD.j^nP J/ + 1 ) = ^ a cm Q 2 cS u 3 I 2 (nD;nP) (2 J' + 1) 

Wji{nP))wj{nD)Sjij , 

where u; is the photon energy, a em is the electromagnetic fine structure constant. In eq.(^) 
one uses 

m c Q 5 - m fe Q c 

Qeff = ; , (6) 

m c + mi, 

where Q Cj b are the electric charges of the quarks. For the B c meson with the parameters from 
the Martin potential, one gets Q e s = 0.41. wj(nL) is the probability that the spin S — 1 in 
the nL state. Sjj> are the statistical factors. The I(nL; nL') value is expressed through the 
radial wave functions, 

I(nL;nL') = | J R nL (r)R nL ,(r)r 3 dr\ . (7) 
For the dipole magnetic M 1-transitions one has 

r(n 1 5 ( - n'Sf + 7 ) = - /i c 2 ff uj 3 (2f + 1) A% , (8) 

where 

A if = J Rns(r)Rns(r)jo(cor/2)r 2 dr , 

and 

A^eff = - (Q c m b - Q- b m c ) . (9) 

2 2m c mb 

Note, in contrast to the ip and T particles, the total width of the B* meson is equal to the 
width of its radiative decay into the B c (0~) state. 

Thus, below the threshold of decay into the BD-pair the theory predicts the existence of 
16 narrow (be) states, which do not annihilate due to the strong interactions, but they have 
the cascade radiative transitions into the ground long-lived pseudoscalar state, the B+ meson. 
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Table 3: The total widths of excited bound (fec)-states below the threshold of decay into the 
-B-D-pair in the model with Martin potential and the branching ratios of the dominant decay 
modes. 



state 


r to t, KeV 


dominant decay 
mode 


BR, % 




0.06 


l^n + T 


100 


2 1 Sn 


67.8 


l^Sn + 7171 


74 


2 1 S 1 


86.3 


+ Try]- 


58 


2 1 Pn 


65.3 


l 1 ^ + 7 


100 


2P 1+ 


89.4 


l 1 ^ + 7 

x *^1 1 / 


87 


2P 1'+ 


139.2 


l 1 5'o + 7 


94 


2 3 P 2 


102.9 


l^i+7 


100 


3^0 


44.8 


2 1 S 1 + 1 


57 


3P 1+ 


65.3 


2 1 S' 1 + 7 


49 


3P 1 /+ 


92.8 


2^0 + 7 


63 


3 3 P 2 


71.6 


2 1 S 1 + 1 


69 


3D 2- 


95.0 


2P 1+ + 7 


47 


3 5 D 3 


107.9 


2 3 P 2 + 7 


71 


3 3 D, 


155.4 


2^0 + 7 


51 


3D 2'- 


122.0 


2P 1 ,+ + 7 


38 



As for the states lying above the threshold of decay into the heavy meson BD pair, the 
width of P* + (3>S) — > B + D°, for example, can be calculated in the framework of sum rules for 



the meson currents, where the scaling relation ]T3[ takes place for the g constants of similar 
decays of T(4S) -> P+P~ and ^(3770) -> D+D~, 



T_ /4mi2\ 
M\M)' 



const. 



The relation is caused by the dependence of energy gap between the vector and pseudoscalar 
heavy meson states: AM 12 • M 12 = const., where Mi >2 are the meson masses in the final 
state, vnyi is their reduced mass. The width of this decay has a strong dependence on the 
B* + (3S) mass, and at M = 7.25 GeV it is equal to T = 90 ± 35 MeV, where the uncertainty 
is determined by the accuracy of method used. 



Table 4: The predictions of scaling relation in comparison with the current experimental data 



value 


exp. 


scaling rel. 


9r(4S)-+B+B- 


52 


input 


<?^(3770)->£>+£>- 


31 


31 


9b*+(3S)^B+D° 




49 
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3. Leptonic constants 



In the framework of potential models the asymptotic behaviour at r — > oo, r — > are 
significant for the determination of the leptonic coupling constants / for the nS'-levels. In 
the leading approximation, the / value does not depend on the spin state of the level and 
it is determined by the value of the radial wave function at the origin, R(0), being model- 
dependent, 



fn 



-Rns{0) 



Taking into account the hard gluon corrections, the constants of vector and pseudoscalar states 
equal 

~ '""- m - (10) 

and the QCD coupling 



fV,P _ f 
J n J n 



1+^ 

7T 



where m 12 are the quark masses, 5 V = 8/3, 5 P = 2 
constant is determined at the scale of the quark masses. 



m 2 



15, 16 



The corresponding uncertainty due to the model dependence is expressed by a factor of 



two. 



Table 5: The radial wave functions at the origin, R n s(Q) (in GeV 3 / 2 ) and R' nP (0) (in GeV 5 / 2 ), 
obtained in the Schrodinger equation with the Martin and BT potentials as well as in the sum 
rules. 



n 


Martin 


BT 


SR 


Ris{0) 


1.31 


1.28 


1.20 


R2S(0) 


0.97 


0.99 


0.85 


R' 2P (0) 


0.55 


0.45 




R' 3P (0) 


0.57 


0.51 





The QCD sum rules |T7[ allow one to determine the leptonic constants for the heavy 
quarkonium states with a much better accuracy. 

Standard schemes of the sum rules give an opportunity to calculate the ground state 
constants for vector and pseudoscalar currents with the account for corrections from the quark- 
gluon condensates, which have the power form over the inverse heavy quark mass. 

QUASILOCAL SUM RULES 

There is a region of the momentum numbers for the spectral density of the two-point 
current correlator, where the condensate contributions are not significant. In this region, the 
integral representation for the contribution of resonances lying bellow the threshold of decay 
into the pair of heavy mesons, allows one to use the regularity of the quarkonium state density 



mentioned above, and to derive the scaling relations [jn| for the leptonic constants of the 
ground state quarkonia with different quark contents and for the excited states. Thus, for 
vector states we have 

fn ( M n\ 2 ( mi+m 2 y _ c (n) 



M n \Mi) V 4m 



1 12 
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where m 12 = m 1 m 2 /(m 1 + m 2 ) is the reduced mass of quarks, and the constant c is determined 
by the average kinetic energy of quarks, the QCD coupling constant at the scale of average 
momentum transfers in the system and the hard gluon correction factor to the vector current. 
Numerically, in the method accuracy the c value turns out to depend on no quark flavour 
and excitation number in the system. Relation ( |TTD is in a good agreement with the data on 
the coupling constants for the families of ip and T particles and, thus, it can be a reliable 
basis for the prediction of leptonic constants for the 5 c -meson family. The constants for the 
pseudoscalar states of nS'-levels are determined by the relation 




1 2 3 4 5 6 7 n 



Figure 2: The dependence of leptonic constants for the nS-levels of bottomonium and char- 
monium, f n , as it is calculated in the Quasilocal Sum Rules. 



Table 6: The leptonic constants for the vector and pseudoscalar states of riS-levels in the (bc)- 
system, f n and ft, calculated in the sum rules resulting in the scaling relation. The accuracy 
is equal to 6%. 



n 


fn, MeV 


ft, MeV 


1 


385 


397 


2 


260 


245 
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4. B c DECAYS 



Lifetime. 

The £? c -meson decay processes can be subdivided into three classes: 

1) the 6-quark decay with the spectator c-quark, 

2) the c-quark decay with the spectator 6-quark and 

3) the annihilation channel — > l + ui(cs,us), where I = e, /x, r. 

In the b — *> ccs decays one separates also the Pauli interference with the c-quark from the 
initial state. In accordance with the given classification, the total width is the sum over the 
partial widths 

T(B C -> X) = T(b -> X) + T(c -> X) + r(ann.) + T(PI) . (12) 

For the annihilation channel the T(ann.) width can be reliably estimated in the framework 
of inclusive approach, where one take the sum of the leptonic and quark decay modes with 
account for the hard gluon corrections to the effective four-quark interaction of weak currents. 
These corrections result in the factor of a% = 1.22 ± 0.04. The width is expressed through the 
leptonic constant of fs c = ff ~ 400 MeV. This estimate of the quark-contribution does not 
depend on a hadronization model, since the large energy release of the order of the meson mass 
takes place. From the following expression one can see that one can neglect the contribution 
by light leptons and quarks, 

r(ann.) = £ ^\V bc \ 2 f B Mm?(l - m?/m| c ) 2 ■ Q , (13) 

i=T,C 

where C T = 1 for the r + ^-channel and C c = 3| V^ s | 2 af for the cs-channel. 

As for the non-annihilation decays, in the approach of the operator product expansion for 
the quark currents of weak decays [|T^], one takes into account the a s -corrections to the free 
quark decays and uses the quark-hadron duality for the final states. Then one considers the 
matrix element for the transition operator over the bound meson state. The latter allows one 
also to take into account effects caused by the motion and virtuality of decaying quark inside 
the meson because of the interaction with the spectator. In this way the b — ► ccs decay mode 
turns to be suppressed almost completely due to the Pauli interference with the charm quark 
from the initial state. Besides, the c-quark decays with the spectator 6-quark is essentially 
suppressed in comparison with the free quark decays because of a large bound energy in the 
initial state. 

In the framework of exclusive approach it is necessary to sum widths of different decay 



modes calculated in the potential models [O, 20]. While considering the semileptonic decays 



due to the b — > cl + vi and c — > sl + ui transitions one finds that in the former decays the hadronic 
final state is practically saturated by the lightest bound 15-state in the (cc)-system, i.e. by 
the r] c and J/ip particles, and in the latter decays the lS-states in the (Es)-system, i.e. B s and 
B*, only can enter the accessible energetic gap. The energy release in the latter transition is 
low in comparison with the meson masses, and, therefore, a visible deviation from the picture 
of quark-hadron duality is possible. Numerical estimates show that the value of B c —>■ (bs)l + vi 
decay width is two times less in the exclusive approach than in the inclusive method, though 
this fact can be caused by the choice of narrow wave package for the B^ mesons in the quark 
model, so that fs 3 ~ 150 MeV, while in the limit of static heavy quark, one should expect 
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a larger value for the leptonic constant. This increase will lead to the widening of the wave 
package and, hence, to the increase of the overlapping integral for the wave functions of B c 
and SW. 

Further, the b — > cud channel, for example, can be calculated through the given decay width 
of b — > cl + vi with account for the color factor and hard gluon corrections to the four-quark 
interaction. It can be also obtained as a sum over the widths of decays with the (uci)-system 
bound states. 

Table 7: The branching ratios of the B c decay modes calculated in the framework of inclusive 
approach and in the exclusive quark model with the parameters |V& C | = 0.040, fs c = 0.47 
GeV, U = 0.54 GeV, f Bs = 0.3 GeV, m b = 4.8 - 4.9 GeV, m c = 1.5 - 1.6 GeV, m s = 0.55 
GeV. The accuracy is about 10%. 



B c decay mode 


Indus., % 


Exclus., % [0 




b — > cl + ui 


3.9 


3.7 


b — > cud 


16.2 


16.7 


Eb^c 


25.0 


25.0 


c — > sl + ui 


8.5 


10.1 


c — > sud 


47.3 


45.4 




64.3 


65.6 


B+ -> T+U T 


2.9 


2.0 


B+ -> cs 


7.2 


7.2 



The results of calculation for the total B c width in the inclusive and exclusive approaches 
give the values consistent with each other, if one takes into account the most significant 
uncertainty related with the choice of quark masses (especially for the charm quark), so that 
finally we have 

r(£+) = 0.55 ±0.15 ps, (14) 

so that the observed B c candidates in the ipn mode at LEP have quite a close value of the 
lifetime. 

Exclusive decays. 

The consideration of exclusive £> c -decay modes supposes an introduction of model for the 
hadronization of quarks into the mesons with the given quantum numbers. The QCD sum 
rules for the three-point correlators of quark currents and potential models are among of those 
hadronization models. 

A feature of the sum rule application to the mesons containing two heavy quarks, is the 
account for the significant role of the coulumb-like a s /t>-corrections due to the gluon exchange 
between the quarks composing the meson and moving with the relative velocity v. So, in the 
semileptonic decays of B£ — ► ^(r/ c )/ + ^, the heavy (Q1Q2) quarkonium is present in the both 
initial and final states, and, therefore, the contribution of coulumb-like corrections exhibits 
in a specially strong form. The use of tree approximation for the perturbative contribution 
into the three-point correlator of quark currents leads to a large deviation between the values 
of transition form-factors, calculated in the sum rules and potential models, respectively [[22 
The account for the a s /t>-corrections removes this contradiction p3| . 
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Thus, the meson potential models, based on the covariant expression for the form-factors 
of weak B c decays through the overlapping of quarkonium wave functions in the initial and 
final states, and the QCD sum rules give the consistent description of semileptonic 5 c -meson 
decays. 

Further, the hadronic decay widths can be obtained on the basis of assumption on the 
factorization of the weak transition between the quarkonia and the hadronization of products 
of the virtual W / * + -boson decay ]2"1 |. The accuracy of factorization has to raise with the 
increase of IV-boson virtuality. This fact is caused by the suppression of interaction in the 
final state. In this way, the hadronic decays can be calculated due to the use of form-factors for 
the semileptonic transitions with the relevant description of W* transition into the hadronic 
state. 



Table 8: The branching ratios of exclusive B c decay modes, calculated in the framework of 
covariant quark model with the parameters |V& C | = 0.040, fs c = 0.47 GeV, /</, = 0.54 GeV, 
f Bs = 0.3 GeV, m b = 4.8 - 4.9 GeV, m c = 1.5 - 1.6 GeV, m s = 0.55 GeV. The accuracy 
equals 10%. 



B c decay mode 


BR, % 


B c decay mode 


BR, % 




2.5 




1.2 


B*l+u t 


6.2 


B a l+v t 


3.9 


lfj7l + 


0.2 


1] c 7l + 


0.2 




5.2 


B S 7T+ 


5.5 


ipp+ 


0.6 


VcP + 


0.5 


B* sP + 


22.9 


B sP + 


11.8 



More generally, the effective four-fermion hamiltonian for the nonleptonic decays of the c- 
and 6-quarks has the form 

K s = ^V uqi V:jCUp)Ol + C c _(p)O c _} + h.c., (15) 

n b cS = ^v qib v^cX{p)oX + c^)ot} + h.c., (16) 

where 

0± = (gia7v(l - 7 5 )c^)('U 7 7 I/ (l - 75)52,5) (<W<V ± S aS 6 j/3 ), 

O b ± = ((710,7,(1 - 75)M (?3 7 7 i/ ( 1 - 75)<?2<5) ((W<5 7 <5 ± S aS S^). 

The factors C± (//) account for the strong corrections to the corresponding four-fermion 
operators because of hard gluons. 

In the leading logarithm approximation at p > m c , one has 

C c^ la^M^)^ 23 fa s (m 2 b ^ ( ' /2 '* 



a s (ml) J V a s (/j, 2 

CQj.) = [c c + (»T 2 . (17) 
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At jj, > rrib, one finds 



Oi s {jn 2 b ) J V a a {n 



T 2 \ x -12/23 / ( lrrt '2\\ -12/25 
-ib ' 



The ai and 02 coefficients, accounting for the renormalization of the four-fermion operators, 
are defined in the following way 

r Nc + 1 +c Nc ~ 1 (20) 

a - C Nc + 1 C Nc ~ l (21) 

In the limit N c — > 00 one has 

01 « 0.5-(C+ + C_), 

a 2 « 0.5 • (C+ — (7_). (22) 

In the theory 02/01 ~ 0.2, which is in agreement with the experimental data, so that the B c 
decay modes proportional to a 2 are essentially suppressed. 




Figure 3: The diagrams of B c 
between the constituents. 



ip(7] c )7T + (p + ) decays with account for the hard gluon exchange 



A decrease of the invariant mass for the hadron system results in an increase of the recoil 
meson momentum. This causes the problem of applicability for the formalism of overlapping 
for the quarkonium wave functions, because, in this kinematics, the narrow wave packages are 
displaced relative to each other in the momentum space into the range of distribution tails. 
In this situation one has to take into account a hard gluon exchange between the quarkonium 
constituents, which destroys the spectator picture of weak transition in the potential approach. 



The widths of 



^7r + and — > 7] c tt + decays E5| have the following forms 



> ljj7T~ 



G 2 F I Vbc 1 2 



128^ /M + 



m 



M 



m 



m) 2 m 



T{B+ -> ^ 



jlfkfW 

(M 

[5(M - m) 2 (M + m) + (M - mf + 8m 3 ] 
16(M + m) 2 M 4 



2 a l 



31 2 



(23) 
(24) 
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where m = or m = m Vc , respectively. 
Numerically one finds 



BR HS (S+ -> 4>7v + ) = 0.77 ±0.19%, (25) 

BR HS (£?+ — > r] c n + ) = 1.00 ±0.25%, (26) 

BR HS (S+ -> i)p + ) = 2.25 ±0.56%, (27) 

BR HS (B+ -> r] c p + ) = 2.78 ±0.70%, (28) 

so that the matrix element is twice enhanced in comparison with that of the potential model 
value due to the contribution of second t-channel diagram. 

The relative yield of excited charmonium states can be also evaluated 

BR"W^V,( 2 S), + )^BR"^^, C (25), + ) 

BR (£?+ — * ipn + ) BR (B+ — > r] c n + ) V ; 

As for the extraction of B c signal in the hadronic background, the decay modes with if) in 
the final state are the most preferable, because the latter particle can be easily identified by 
its leptonic decay mode. This advantage is absent in the B c decay modes with the final state 
containing the r/ c or B^ mesons, which are the objects, whose experimental registration is 
impeded by a large hadron background. 

From the values of branching ratios shown above one can easily obtain, that the total 
probability of ip yield in the B c decays equals BR(B+ — > ipX) = 0.17. 

It is worth to note that the key role in the B c signal observation plays the presence of the 
vertex detector, which allows one to extract events with the weak decays of long-lived particles 
containing heavy quarks. In the case under consideration it gives the possibility to suppress 
the background from the direct ip production. In the semileptonic — > ipl + vi decays the 
presence of vertex detector and large statistics of events allows one to determine the S c -meson 
mass and to separate the events with its decays from the ordinary i?+-meson decays, which 
have no mode. In the if)7i + decay the direct measurement of B c mass is possible. The 
detector efficiency in the reconstruction of three-particle secondary vertex from decays 

of ip and 7T + or /+) becomes the most important characteristics here. The low efficiency of the 
LEP detectors (e as 0.15), for example, makes the B c observation to be hardly reachable in 
the experiments at the electron-positron collider. 

5. B c PRODUCTION 



The (be) system is a heavy quarkonium, i.e. it contains two heavy quarks. This determines 
the general features for the B c meson production in various interactions: 

1. Perturbative calculations for the hard associative production of two heavy pairs of cc 
and bb cause the suppressed B c yield as a 10 -3 fraction of beauty hadrons. 

2. A soft nonperturbative binding of nonrelativistic quarks in the color-singlet state can 
be described in the framework of potential models. 

The latter means the factorization of hard and soft amplitudes, so that the soft factor is 
determined by the radial wave function at the origin for the S-wave quarkonium and its first 
derivative for the P-wave one. Thus, the analysis of B c production mechanisms is reduced to 
the consideration of perturbative amplitudes in high orders over a s . 



14 



e + e annihilation 



In e + e~-annihilation at large energies (M§ /s <C 1), the consideration of leading diagrams 
for the 5 c -meson production gives the factorized scaling result for the differential cross-section 
over the energy fraction carried out by meson, da/dz = a(bb) ■ D(z), where z = 2Es c /y/s. 
This distribution allows the interpretation as the hard production of heavier 6-quark with 
the subsequent fragmentation into B c , so that D(z) is the fragmentation function. In this 
approach one can obtain analytic expressions for the fragmentation functions into the different 
spin states of S, P, and .D-wave levels Eft . 




B r 






B,. 



Figure 4: The diagrams of the 5 c -meson production in e + e -annihilation 
Denoting r = m c j [m^ + m c ), for the pseudoscalar state one finds 



8a 2 \^(0)\ 2 rz(l-z) 2 . 9N 9 

A 1 (6 - 18(1 - 2r)z + (21 - 74r + 68r 2 )z 2 - 



Mml (l-(l-r)z) 6 
2(1 - r)(6 - 19r + 18r 2 )z 3 + 3(1 - r) 2 (l - 2r + 2r 2 )z 4 ) 

for the vector meson the perturbative fragmentation function is equal to 

8a 2 |^(0)| 2 rz(l-z) 2 



(30) 



d ^ = 1 (i - (i - ;w (2 - 2(3 - 2r)z + 3(3 - 2r + 4r v 

2(1 - r)(4 - r + 2r 2 )/ + (1 - r) 2 (3 - 2r + 2r 2 )z 4 ), 



(31) 
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These expressions take into account the spin structure of both interactions and bound 
states, but they are very close to the Peterson et al. model |27|, where 



1 1 rz(l — z 

D i Z h-,B^ ~ -7— 2 Ap 



^K-^-S) 2 (l-(l-r)^) 4 ' 

so that the calculated functions are slightly more hard. Numerically, one finds f(b—> B£) = 
c(B+) I cr(bb) = (1.3 ± 0.4) ■ 10 -3 . Combining the latter with the prediction of branching ratios 
one gets 

[f{b -> B+) ■ BR(B+ -> ^vr + )] TH = (0.22 ± 0.09) • 10~ 5 , 

where f(b — > 5+) = o~(Z — > B^)/a(Z — > gg), which may be compared with the OPAL 
estimate 

[f{b - 5+) • BR{B+ - ^ + )]opal = (3.8±i:2 ± 0.5) ■ lO" 5 . 

If one takes into account both the single DELPHI candidate in the ipir mode and absence of 
the event at ALEPH, then assuming the equal detector-efficiencies, one estimates 

[/(& ^ B +) ■ BR{B+ - ^7T+)] LEP = (1.9lf:l ± 0.3) • 10- 5 . 

Photon-photon collisions 

In photon-photon interactions for the B c production in the leading approximation of per- 
turbation theory, one can isolate three gauge-invariant groups of diagrams, which can be 
interpreted as 

1. the hard photon-photon production of bb with the subsequent fragmentation of b — > 
B£(nL), where n is the principal quantum number of the (6c)-quarkonium, L — 0, 1 ... is the 
orbital angular momentum, 

2. the corresponding production and fragmentation for the c-quarks, and 

3. the recombination diagrams of (6c)-pair into B+, wherein the quarks of different flavours 
are connected to the different photon lines. 

In this case, the results of calculation for the complete set of diagrams in the leading order 
of perturbation theory show that the group of ^-fragmentation diagrams at high transverse 
momenta pt{B c ) ^> M# c can be described by the fragmentation model with the fragmentation 
function D^ B +(z), calculated in the e + e~-annihilation. The set of c- fragmentation diagrams 
does not allow the description in the framework of fragmentation model. The recombination 
diagrams give the dominant contribution to the total cross-section for the photon-photon 



production of B c [28 



The reason is that the photon coupled to the charm quark can be considered as the "re- 
solved" one, so that the bottom quark is hardly produced off the charm quark appearing in 
the structure function of the photon, as it is evaluated in the lowest order of perturbative 
theory. 



16 



da e + e -/dz, pb 

0.18 
0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 5: The differential cross-sections over z for the B c production in e + e~ — > 7* — > B c + 
X (histograms) in comparison with the fragmentation contribution calculated analytically 
(smooth curves) at 100 GeV for the following states of B c : 1 P\ (solid line), 3 Pq (dashed line), 
3 Pi (dotted line), 3 P2 (dash-dotted line). 
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Figure 6: The contribution of fragmentation for b — > 5+ in the photonic production (three 
initial diagrams) and the recombination (the last diagram). 
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Figure 7: The summed total cross-section versus the energy of photonic production for the 
P-wave states of B c (•) in comparison with the fragmentation contribution (dashed line). The 
solid-line curve is the fit of points. 
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da-yy/dz, fb da.yy/dz, fb 
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' 3 p ' ' 3 p 

-r 1 Z *2Z 

Figure 8: The differential distributions over z for the photonic production of P-wave 
states of B c at 100 GeV: complete set of diagrams (solid histogram), the 6-fragmentation 
diagrams in comparison with the prediction of fragmentation model (dashed histogram 
and curve), the c-fragmentation diagrams and the fragmentation model (dotted histogram 
and curve), correspondingly. 
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Figure 9: The distributions over the transverse momentum. 
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Hadronic production of B c 



The parton subprocess of gluon-gluon fusion gg — > B+ + b + c dominates in the hadron- 
hadron production of B c mesons. In the leading approximation of QCD perturbation theory 
it requires the calculation of 36 diagrams in the fourth order over the a s coupling constant. 
In this case, there are no isolated gauge-invariant groups of diagrams, which would allow the 
interpretation similar to the consideration of B c production in e + e~-annihilation and photon- 
photon collisions. 

By the general theorem on factorization, it is clear that at high transverse momenta the 
fragmentation of the heavier quark Q — > (Qq) + q must dominate. It is described by the 
factorized formula 

da_ = l d»( K9g -.QQ) D x . d, 

dp T J akx x 

where [i is the factorization scale, da/dkx is the cross-section for the gluon-gluon production 
of quarks Q + Q, D is the fragmentation function. 

The calculation for the complete set of diagrams of the 0(af)-contribution allows one 
to determine a value of the transverse momentum pjt m , being the low boundary of the region, 
where the subprocess of gluon-gluon 5 c -meson production enters the regime of factorization for 
the hard production of 66-pair and the subsequent fragmentation of 6-quark into the bound 
(frc)-state, as it follows from the theorem on the factorization of the hard processes in the 
perturbative QCD. 

The p™ m value turns out to be much greater than the M# c mass, so that the dominant con- 
tribution into the total cross-section of gluon-gluon B c -pro duct ion is given by the diagrams of 
nonfragmentational type, i.e. by the recombination of heavy quarks. Furthermore, the convo- 
lution of the parton cross-section with the gluon distributions inside the initial hadrons leads to 
the suppression of contributions at large transverse momenta as well as the subprocesses with 
large energy in the system of parton mass centre, so that the main contribution into the total 
cross-section of hadronic _B c -pro duct ion is given by the region of energies less or comparable 
to the S c -meson mass, where the fragmentation model can not be applied by its construction. 
Therefore, one must perform the calculations with the account for all contributions in the 
given order under consideration in the region near the threshold. 

The large numeric value of p™ in points out the fact that the basic amount of events of 
the hadronic .B^-production does not certainly allow the description in the framework of 
the fragmentation model. This conclusion looks more evident, if one considers the B c - meson 
spectrum over the energy. 

The basic part of events for the gluon-gluon production of B c is accumulated in the region 
of low z close to 0, where the recombination being essentially greater than the fragmentation, 
dominates. One can draw the conclusion on the essential destructive interference in the region 
of z close to 1, for the pseudoscalar state. 

We have in detail considered the contributions of each diagram in the region of z — > 1. In 
the covariant Feynman gauge the diagrams of the gluon-gluon production of Q + Q with the 
subsequent Q — > (Qq) fragmentation dominate as well as the diagrams, when the qq pair is 
produced in the region of the initial gluon splitting. However, the contribution of the latter 
diagrams leads to the destructive interference with the fragmentation amplitude, and this 
results in the "reduction" of the production cross-section in the region of z close to 1. In the 
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Figure 10: The differential cross-section for the meson production in gluon-gluon collisions 
as calculated in the perturbative QCD over the complete set of diagrams in the 0(af) order 
at 200 GeV. The dashed and solid histograms present the pseudoscalar and vector states, re- 
spectively, in comparison with the results of fragmentation model shown by the corresponding 
smooth curves. 
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Figure 11: da/dz for the B c - mesons at the gluon interaction energy of 100 GeV. The dashed 
histogram presents the gg — > B c + b + c process, the dotted one is the abelian case. The 
curve shows the result of the fragmentation model. The cross-sections are normalized over the 
cross-section of the bb pair production. 
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datot/dpT, pb/GeV 




GeV 

Figure 12: The differential cross-sections for the B c (B*) -meson production in pp-collisions 
versus the transverse momentum at the interaction energy 1.8 TeV with the cut off the gluon- 
gluon energy E > 60 GeV. 

axial gauge with the vector = Pq this effect of the interference still manifests itself brighter, 
since the diagrams like the splitting of gluons dominate by several orders of magnitude over the 
fragmentation, but the destructive interference results in the cancellation of such extremely 
large contributions. This interference is caused by the nonabelian nature of QCD, i.e., by the 
presence of the gluon self-action vertices. 

To stress the role of the interference diagrams related to the nonabelian self-action of 
gluons, we have considered the process with abelian currents. In the abelian case the effect 
of the destructive interference due to the additional contribution of the self-action of gauge 
quanta, is absent. So, the agreement between the factorized model of fragmentation and the 
exact perturbative calculation is quite good at z close to 1. 

The direct verification of the given mechanism for the Y3 c -meson production could be the 
comparison of the 73 c -meson spectra in two hemispheres in the region of the gluon fragmenta- 
tion and in the photon one, in the photonic production of B c on nucleons. 



At LHC with the luminosity C = 10 34 cm _2 s _1 and y/s = 14 TeV one could expect 

4.5 ■ 10 10 B^ events per year, 
that leads at the efficiency of the decay reconstruction e = 0.1 to 

3 • 10 4 B^ — > ^7r + events per year, 
with no taking into account some detector-design cuts off angles and momenta. 
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6. Conclusions 



The family of (be) mesons contains 16 narrow states, the excited ones decay into the 
ground pseudoscalar state due to the radiative cascades. 

The mass of ground state is expected to have the value 

m Bc = 6.25 ± 0.03 GeV, 
which is quite close to the value of B c — > ijm candidates at LEP 

KbJlep = 6.33 ± 0.05 GeV. 

The meson is the long-lived particle with the predicted lifetime equal to 

r Bc = 0.55 ±0.15 ps, 
which must be compared with the LEP measurements of the candidates 

[tbJlep = 0.281H ps. 

In e + e~-annihilation, the fragmentation of b — > dominates, so that one expects 

[f(b -> B+) ■ BR(B+ -> ^7t + )] th = (0.22 ± 0.09) • 10" 5 , 
which is less than the B+ — > ^7r + candidates rate at LEP 

[f(b ^ B+) ■ BR(B+ ^ ^7T + )] LEP = (1.9±?;| ± 0.3) • 10" 5 . 

The fragmentation regime works in the hadronic production of B c at high transverse 
momenta pr > 35 GeV only. The recombination is essential at lower momenta. In 
the forward production of B c along the hadron beam, the destructive interference of 
fragmentation diagrams with the gluon-splitting ones takes place. 

The most promising production rate of B c for its observation is expected in hadron- 
hadron collisions, especially at LHC and upgraded Tevatron. 
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